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. In this review, I describe a recent approach based on the representation theory of the Onsager algebra which 

^ I is used to derive exact resuhs for the XXZ open spin chain. The complete spectrum and eigenstates are obtained 
as rational functions of a single variable which discrete values correspond to the roots of a certain characteristic 

' ^ ' polynomial. Comments and open problems are also presented. 
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■ 1 Introduction 

I ■ 

■ Since Sklyanin's seminal work [T] on spin chains with integrable boundary conditions, finding exact results such 
, as the energy spectrum of a model, corresponding eigenstates and correlation functions for elementary excitations 

xl ■ has remained an interesting problem in connection with condensed matter or high-energy physics. Among the 
known integrable open spin chains that have been considered in details, one finds the XXZ open spin chain. For 



\. generic boundary conditions, its Hamiltonian reads 

! N-l 

. k=l 

><■ (,1/2 _,-l/2) (,W ^ 2 
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where cri,2,3 and a± = (ai ± icr2)/2 are usual Pauh matrices and A (q^^^ + q^^^^)/2 denotes the anisotropy 
parameter. Here, we restrict our attention to the (massless) regime — 1 < A < 1 i.e. q = exp(<^) with (p purely 
imaginary. To exhibit the six independent boundary parameters, below we will sometimes use the following 
parametrization with 6,9 eM. and a,a d C: 

= (eL"V = cosha , fc+ = (fc_)t = -(gi/2 _ q-y^)e'0 (left) , 

K* = (k^ ^ - cosh5 , K+ = = -e'V(2(gi/2 + ,-i/2)) ^^^g^it) . (2) 

Known to be integrable since 1988, for the simplest diagonal boundary conditions i.e. k± = k± = the 
algebraic Bethe ansatz approach has been succesfully applied and used to derive the spectrum, eigenstates [1] and 
more recently correlations functions [5] of the model. The derivation of these results essentially relies on the fact 
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that the pseudo- vacuum (or reference) state of this particular case of open spin chain actually coincides with the 
one associated with the XXZ spin chain with periodic boundary conditions. 

For non-diagonal boundary conditions, there is no simple or obvious pseudo- vacuum state. From 2001 to 2007, 
several attempts [SI IH [5j |6] were considered in order to circumvent this difficulty: either based on the functional 
Bethe ansatz approach, or based on the algebraic Bethe ansatz approach using suitable gauge transformations of 
the double-row monodromy matrix provided the boundary parameters satisfy certain linear relations or 5 is a root 
of unity. In any case, either the spectrum was obtained but not the eigenstates or the spectrum and eigenstates 
were obtained but only for certain sets of boundary conditions. The main advantage of a Bethe- type solution being 
its efficiency in studying the thermodynamic limit of the spin chain through the link between Bethe equations and 
non-linear integral equations, it may explain why Bethe ansatz approaches (functional or algebraic) have deserved 
such attention in all attempts at solving the XXZ open spin chain ([T]). However, for generic boundary parameters 
and q, these standard approaches have failed. Note however a rather simple functional approach recently proposed 
for generic parameters, where the spectrum is written in terms of solutions of a system of highly transcendental 
equations which are not Bethe equations 7J . 

Motivated by the challenging problem of solving the XXZ open spin chain ([!]) in the most general regime of 
integrable boundary conditions, in 2005 we decided to explore a different and non-standard (compared to Bethe 
ansatz) path which takes its roots in Onsager's original article on the solution of the planar Ising model [8]. A 
novel approach followed, essentially based on the analysis of the representation theory of the 5— Onsager algebra. 
After some efforts [9l [101 HI] i we used this approach to derive the spectrum and the eigenstates of the XXZ open 
spin chain ([l} for generic boundary conditions and q, as well as we recovered the known results for special relations 
between left and right boundary parameters [12] . 

In this review, our main objective will be to present the novel approach and methodology we have developped 
in the last two years and describe the main results rather than the technical aspects or details which can be 
found in our papers [9l I10 [ fTT | I12j written in collaboration with K. Koizumi. In addition, comments and further 
directions are briefly presented. 



2 The strategy 

The starting point of the analysis was the discovery that the transfer matrix txxz{u) generating ^ can be written 
in terms of a finite set of mutually commuting quantities X21+1 which form a deformation of the Onsager's or 
Dolan-Grady hierarchy [13l[14]. Namely, one has: 

N~l 

txxziu) - ^2i+i{u) I21+1 + Mu) I with J21+1 - kW''^i^ + n*w['^l + ^§["11 + ^gl^^} (3) 

where J^2i+i{u)^ 3^o{u) are rational functions that have been explicitely calculated in (lOJ. In this description, note 
that the remaining boundary parameters in Q enter in the explicit realizations of the AN fundamental operators 
[To] . Such remarkable structure comes from the fact that the elements of the double-row monodromy matrix (also 
called Sklyanin's operator) T(u) - u being the spectral parameter - associated with the XXZ open spin chain ([T]) 
can be decomposed on a basis of AN operators wi^'' , W^^J , , ^[.^J , k = 0, l,...,iV — 1, formally as (details 
can be found in [l0|): 



/ ^(^;{WW,WW}) B{u-{gil\}) I 



In addition, using the fact that the elements of the double-row monodromy matrix satisfy the reflection algebra 
it is possible to show that these (spectral parameter independent) operators - which act on TV-tensor product 
representations of C/„i/2(s?2) - induce a finite dimensional representation of the g— deformed analogue of the 
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Onsager infinite dimensional algebra |9j with defining relations!! for k,l 

[Wo, W,+i] = [W_fe, Wi] = ^-^^^-l_^(G,+i - G,+0 , 

[Wo,Gfc+i]^= [Gfe+i,Wo]^ = pW_fe_i-pWfe+i , 
[Gk+i, Wi] ^ = [Wi , Gfc+i] ^ = pWfe+2 - pW_fe , 
[Wo,W_fc]-0, [Wi,Wfc+i]=0, 

[Gfc+i,G,+i] = , [Gfe+i,G;+i] =0 , [Gfe+i,G,+i] + [Gfe+i,G,+i] ==0 . (5) 

In light of the above formulation, it was tempting to consider the XXZ open spin chain ([1]) from the point of 
view of the representation theory of the algebra ([5]) in view of the success of the Onsager's algebra and related 
representation theory in solving the Ising model |8j, generalizations |14j and superintegrable chiral Potts model 
[151 [TBj . Having this in mind, we implemented the following procedure in order to solve the spectral problem for 
([T]) or, more generally, the transfer matrix: 

• Use the relation between the algebra ([5]) and Terwilliger's tridiagonal algebra [T7]. See [S]; 

• Show that W^'^\}V['^^ form a tridiagonal pair (see below). See P^igifTT]: 

• Construct the states {"^i^]} ^^^d {<y3^j^|} on which Wq^\w[^'' act as block tridiagonal matrices. See [TT] : 

• Show that yi^'^\y\^i^i,Gk+i,Gk+l, k = 1, iV — 1 also act as block tridiagonal matrices; See [H]. 

• Write the spectral problem for all quantities 22/+1 as: 

— > a coupled system three-term recursion relations; 
— > a system of second-order q— difference equations; 
and derive the spectrum and eigenstates. See [T^ . 

It is important to stress that the above procedure introduced to derive all eigenvalues and eigenstates of ^ 
solely uses the block tridiagonal structure of the commuting operators 22/+ 1 in two suitable basis. However, simi- 
larly to the Ising or superintegrable model [TQ , it is possible to show that the operators wi^^ , VV^^j , G[!^\ , gI!^\ , 
fc = 0, ...,N satisfy a set of linear relations. Exploiting this property should lead to another solution - may be 
simpler - to the spectral problem for the transfer matrix of the XXZ open spin chain ([1]). We will explain this 
point in the last Section. 

3 Solving the spectral problem 
3.1 Mathematical background 

Almost nothing being known about the representation theory of the algebra (O, we used its connection with 
Terwilliger's tridiagonal algebra [171 (sometimes called the g— Onsager algebra). Indeed, we had previously noticed 
[TD] (see also [H]) that >V^^\ w]^ given by 

W^^' - (A:+a+ + A:_a_)®j(^-i)+g-^/2®wr"'\ 

Wi""^ = (fc+a+ + A;_a_)® 1(^-1) +g---^/2®>vl'^-i) (6) 
with Wg'''' = e^"* , w[^'' = satisfy the so-called g— Dolan-Grady relations 



^The g-commutator [X,Y]q = q^/^XY - q-^/^YX is introduced and the scalar p = (q^/^ + q-'^^^)^k+k- for the XXZ open 
spin chain. 
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and induce a representation of the tridiagonal algebra. According to Terwilliger's work [T7], it follows that 
Wq^\ w}^' provide an explicit example of tridiagonal pair. Roughly speaking, for generic boundary parameters 
^ and generic values of q it means that there exists a complete basis (resp. 'dual' basis) of V in which Wq^'' 
(resp. w}^'') is represented by a diagonal matrix with degeneracies and w[^^ (resp. Wp^'') is represented by a 
block tridiagonal matrix. Such basis and dual one have been constructed in jTT], and are denoted below {"^^^j} 

and {'/'ijl]}- 

Let us recall the main results of [11 . Using the canonical basis ^'jLi on which the nonlocal operators 
^ act, let us introduce the 2^ ordered states 

<1 - {e--^^-'-'-^'^'-^'\+U + \-U^ for ,e{(^;^)+i,...,(^)}.(8) 

For generic boundary parameters and q, these states form a complete[f| basis of V on which the tridiagonal pair 
Wq^\ W^''^^ acts [11. Using the explicit expressions ([6]) it is straightforward to exhibit the diagonal and block 
tridiagonal structure of the operators 'Wq^\ 'W[^\ respectively, in this basis: 

wfViS = ^^rr¥^fi ^ith AW=cosh(a + (iV-2n)^/2), 

i=l i—1 i—1 

with n = 0, 1, iV and jG{l,...,(^)}. Note that the explicit form of the coefficients ol"^^^^\ ^n[ij]' '^^n 

be found in [11] . Clearly, a similar analysis can also be performed by considering the eigenbasis of w}^' instead 
of Wq^\ leading to analogous results. In this case, the 'dual' states given by (/s^j^j = ■f/'ljlj* U^-q*, (j}^-<t>. e^e+i-n 
form an eig enbasis of W^^^ [TT]. 

All higher operators realizing ([S]) being generated recursively from the lower ones |10| , there were good reasons 
to expect that they would act in a simple manner on the eigenbasis of Wq^' or w[^\ Indeed, proceeding similarly 
for higher values of Z — 1,2, ...,iV — 1 it is possible to show that all nonlocal operators wi^'' , W^^j , , 
also enjoy a block tridiagonal structure in the basis ([5]). After some straightforward calculations, one finds: 

2 = 1 

^^l+l^n[j] - "nfo] 'Pn+lfi] +Z^ "„[y] V„[,] + Z^ '^nfo] V^n-l[i] ' 

i—1 i—1 i—1 

(" + i) (") 
- (n^/^\W o-i/2\W ^7,(^''\/,W , 

i=l i=l 

+ 2^ \1 1 ^n-ll^n[%i]yn-l[i\ ^ 

i=l 

i=l i=l 



^For special relations among the boundary parameters, the representation may become indecomposable. We do not consider such 
possibility here. 
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+ eV/^aS - <f'''>^')C^i^,^ (10) 

where all coefficients can be found in [T^] . Note that they have a rather simple form regarding to the recursion on 
the index I, a fact which may have been anticipated from the existence of linear relations among the operators. 
See the comments below. In the dual basis similar results are also obtained. 

Having two different basis {?/^^^.|}, {(/Sgj^j} on which all mutually commuting operators generating the Dolan- 
Grady hierarchy Q act as block tridiagonal matrices with entries known explicitely, we are now ready to consider 
the spectral problem for the XXZ open spin chain ^ formulated using 

3.2 Eigenstates and spectrum for generic parameters 

For generic values of the boundary parameters and q, there is no nontrivial subspace of V which is left invariant 
under the action of the nonlocal commuting operators T21+1. Considering the most general linear combinations of 
states V'i^] or '/'l|fej; it follows that any eigenstate admits two dual representations with respect to the two dual 
basis on which the elements of the g— Onsager algebra act in a block tridiagonal form. Any eigenstate can be 
written cither 

N {I) N (^) 

n=0 j=l s=0 k=l 

where each of the dual families of weights {/^^j}, {/fjfcj} and the restricted set of 'allowed' eigenvalues {Ai^i, A^ 2"} 
have to be determined. For generic boundary parameters, the spectrum of Ii is not degenerate. As [1^21+1, ^2k+i] = 
0, it follows that {f^j\}^ {/i[fe]} ^e uniquely derived from the constraint: 

XilvI'(Ai)) = Ai|vI/(Ai)) . (12) 

For any representation ()lip choosen, this equation leads to a system of coupled three-term recurrence relations for 
the weights which, in some sense, generalizes the Askey- Wilson recursion relations0 for g— orthogonal polynomials. 
For n, s — 0,1, N one has: 

E ^i-fbm]/«-lH(Al) + E ^l+fbrnl-^n+lH^^O + E ("^ibm] ~ ^l^J™) ^^l) = ^ ' (1^) 
m— 1 m— 1 m— 1 

or, alternatively, 

E^ «I"[l,/i"iri (Ai) + E^ clTiilj/l+ii, (Ai) + E (X£T - A,4™)/itr' (Ai) = (14) 

m— 1 m— 1 m— 1 

where the coefficients '^1^'^^], '^s[kni] ('^ ^ i^' ^^^}) '^^^ easily obtained from (fTU|) [T^ . Having 2^ equations 
in total, one finds that these weights are rational functions of the variable Ai G {Ai_i, A^ 2«} where Ai^^. are 
the roots of the characteristic polynomial of degree d = 2^ 

V{Ki) =dct[Xi - AiJ] . (15) 

The complete family of eigenstates (jlip being identified, the spectrum of the transfer matrix can now be derived. 
For instance, acting with higher operators X2;+i of the g— Dolan-Grady hierarchy on these eigenstates, and taking 
the scalar product with 

^iN) _^(iV), TW _ (TV), 



For N = 1, the coefBcients are special cases of the Askey- Wilson ones. 
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respectively, according to the representation choosen in (jlip one ends up with a second-order q— difference equation 
that determines all higher necessary eigenvalues A21+1, I = 1, N — 1 as rational functions of Ai. For any choice 
oi s,k,n,j^ one writes either 



m=l y^[fc] m=l ISJ m=l [fc] l^j 

m=l ^[j] \^) m=l ^[j] l*^; m=l ^ [j] l^^J 



where the coefficients can be found in [T^ and the functions 

have been introduced. To resume, all eigenvalues I = 1,...,-/V — 1 are obtained as rational functions of 

the single variable Ai restricted on the discrete support defined by P(Ai) — with In terms of these, the 

spectrum of the Hamiltonian ([T]) reads: 



2(gi/2 + g-i/2)^ 

3.3 Truncation and the Bethe ansatz regime of boundary parameters 

For certain special relations between left and right boundary parameters, important simplifications occur. These 
situations arise when some of the off-diagonal lower (resp. upper) blocks of the block tridiagonal matrices repre- 
senting X21+1 for any / = 0, l,...,iV — 1 identically vanish, leading respectively to two different conditions on the 
boundary parameters. For all I, i, j, to, k and integer p fixed: 

bI^^=0 and S;!"^;^,;,,^,] = ^ a±& = -i{0 - 9) - {N - 2p - 1)^/2 mod{2i7r) ; (18) 
CiM^O and c'Zil,,^^0 ^ a±a = i{e - 0) - {N - 2p - 1^2 mod{2in) . (19) 

In these two special cases, the coupled three-term recursion relations (|13p or (|14p determining the weights in the 
expansions (llip become truncated. This means that any eigenstate can be written as a truncated combination of 
one of the two representations in pTjl . For the first family of boundary conditions (US]), the eigenstates read: 

P (") N-p-l{l) 
n=0 j=l s=0 k=l 

whereas 

*^"nAO= E E/Sr'lAO'Aa and vI.(-)(A,)= E E /S^HaO^IS (21) 

n=p+l j = l s=N—p k=l 

for the second family of boundary conditions (|19p . In others words, the complete family of eigenvalues and 
eigenstates splits in two sets, each associated with a characteristic polynomial of degree c?i < 2^ and ^2 < 2^, 
respectively, such that o?i + c?2 = 2^. According to the family of boundary conditions one finds: 



di = E I J for (UHl) and = E ( ) ^"'^ ^ 
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Then, the eigenvalues A21+1 for I — 1, N are obtained with the same method as for generic boundary parameters. 

It is important to mention that the hnear relations (|18p , p9p among the boundary parameters coincide exactly 
with the ones for which a Bethe-type solution exists, although the derivation of these relations is here totally 
different. Also, the splitting of the eigenvalues and eigenstates in two sectors is also confirmed: in the algebraic 
Bethe ansatz framework [H [5] (see also p]), the two sets of eigenstates and eigenvalues are described by two 
different sets of Bethe equations. This shows that the approach described here provides an alternative solution 
compared to the Bethe ansatz one for the model ([T]) with the family of boundary conditions (fTSl) . For these 

regimes of boundary parameters, a detailed comparison between both approaches is an interesting problem: it 
should give explicit relations between the solutions of Bethe equations and the roots of the truncated characteristic 
polynomials. 



4 Comments and open problems 

• Is the model 'rea/Zw' solved? If one believes that solving an integrable model only means that one must 
be able to find its spectruml^in terms of solutions of highly transcendental equations (Bethe-type equations), then 
the solution we proposed for the XXZ open spin chain ([T]) (12j is indeed not satisfying at all. And from this point 
of view, this is not the only integrable model which is out of rangeHl On the other hand, if one thinks that solving 
an integrable model means that one is able to reduce the complexity of the system using unified mathematical 
structures (quantum algebra, representation theory) in order to write its spectrum in terms of solutions of simpler 
equations (for instance the algebraic equations with (I13|) . (|14|) where Ii is a block tridiagonal matrix) then from 
this point of view the model ^ is solved. 

• Large N and thermodynamic limit? For generic boundary parameters, a natural question to ask is 
weither diagonalizing Xi is simpler and faster than diagonalizing TL and all higher conserved quantities generated 
from the transfer matrix. Indeed, one has to remember that the eigenstates and eigenvalues of Ii generate the 
eigenvalues of all higher quantities 22/+ 1 for Z = l,...,iV— 1. Thanks to the block tridiagonal structure of the 
matrix Zi, finding the spectrum has been reduced to (fT3|) or (fT4|) . Although this simplified formulation helps for 
numerical calculations (using well-known FORTRAN procedures for block tridiagonal matrices), for iV >> 1 or in 
the thermodynamic limit N —* 00 solving directly these 2^ coupled three-term recurrence relations is a problem 
that clearly needs further investigations. Studying infinite dimensional representations of ([5]) might be a possible 
direction. 

• Another solution using the Onsager algebra? The solution here described [T2] is essentially based 
on the block tridiagonal structure of all operators satisfying ([5]) in a certain basis. However, for finite dimensional 
representations associated with the XXZ open spin chain ^ it is also interesting to point out that0 

1=1 

N 

1=1 



(qi/2 _ q- 


-1/2) 


k-\-k- 




(5I/2 _ q- 


-1/2) 


k-^-k^ 






-1/2) 


k-f-k^ 




(q'/'-q- 


-1/2) 



1=1 
N 



rgr + Y.ci'!liOi:i = o . (22) 



k+k^ 

+ 1=1 



®I do not even speak about the eigenstates which can not be derived using the functional Bethe ansatz approach. 
^Consider for instance the XYZ chain, quantum Toda chain, chiral Potts model,... 

/ , N-n-2 / . (0) (0), 1/2 -1/2n2 

^-n -\ '-J Vy -l-y ) \^(ql/2+g-l/2) J (u+1) ! {iV-n-2) ! I (ql/2+q-l/2) JV_„_1 "t" fe^ fc_ (gl/2 + g- 1/2 ) 
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For the reader who is famihar with the Ising and superintegrable chiral Potts models, the spectrum of these 
models is known to be expressed in terms of the roots of the polynomial (called Baxter's polynomial for the 
superintegrable Potts model) associated with '^iCtiAi = and aiGi = where Ai,Gi are Onsager's algebra 
generators (see references and details in [E]). By analogy, linear relations of the form ([22|) may be key ingredients 
in finding another solution - may be simpler than the one we proposed in [12| - to the XXZ open spin chain ([1]). 

• Beyond. The computation of correlation functions essentially relies on the knowledge of the complete set of 
eigenstates (in the massless regime which is here considered). In view of (jlip . one can easily calculate the action 
of the elementary local operators on the basis ([8]) 

N iZ) 

^i'^^S-EE-^rji^M^S ^it^ ^ a e {1,2, 3}, 

and the coefficients A^^^^j^j ^j^j are derived recursively. As an example, it is possible to write one-point functions 
as rational functions of the variable Ai defined on the discret support V{Ai) = with (fT5|) . For small values 
of A'' < 6 and generic parameters or the (Bethe ansatz) regime of parameters!! (|18|) . (fT9|) . the results are easily 
checked numerically. Having a better understanding of solutions to (fT3|) . (fT4|) at large A^ would then play an 
essential role for studying correlation functions. In this direction, extending the results of lTDJ might be helpful. 

Acknov^rledgements: I thank the organizers of RAQIS Meeting 2007 where these results have been presented. 
Part of this work is supported by the ANR research project ^''Boundary integrable models: algebraic structures 
and correlation functions" , contract number JC05-52749. 
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